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Abstract In recent years, a number of experimental studies have been conducted 
to investigate the mechanical behavior of water-saturated biological tissues like 
articular cartilage under dynamic loading. For in vivo measurements of tissue via- 
bility, the indentation tests with the half-sinusoidal loading history were proposed. 
In the present paper, the sinusoidally-driven compression test utilizing either the 
load-controlled or displacement-controlled loading protocol are considered in the 
framework of linear biphasic layer model. Closed-form analytical solutions for the 
integral characteristics of the compression test are obtained. 
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1 Introduction 



It is well known that the long-term creep and relaxation tests, typically used for 
determining viscoelastic and biphasic/poroelastic properties, are not appropriate 
for rapidly assessing the dynamic biomechanical properties of biological tissues 
like articular cartilage. For in vivo measurements of tissue viability, [Appleyard et| 



al. (20011 developed a dynamic indentation instrument, which employs a single- 



frequency (20 Hz) sinusoidal oscillatory waveform superimposed on a carrier load. 
It is to note that the oscillation test requires some time period to elapse before 
recording the measurement to minimize the influence of the initial conditions. That 
is why the indentation or compression tests with the half-sinusoidal loading history 
are so tempting. In particular, the half-sinusoidal history is useful for developing 



indentation- type scanning probe-based nanodevices reported by Stolz et al. (20071. 



Also, as a first approximation, such a indentation history can be used for modeling 



impact tests (Butcher and Segalman[ 1984). 
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In the present paper, it is assumed that the mechanical response of a time- 
dependent material can be described in the framework of biphasic layer model. 



Following Argatov (20121, we consider sinusoidally-driven flat-ended compression 



test utilizing either the load-controlled or displacement-controlled loading protocol. 
Closed-form analytical solutions for the integral characteristics of the compression 
test are obtained. 



2 Unconfined compression of a cylindrical biphasic sample 



We assume that the unconfined compression test for a soft biological tissue sam- 
ple can be described by a one-dimensional axisymmetric mathematical model of 



Armstrong et al. (1984) developed in the framework of the linear biphasic theory 
(Mow et al. 1980). In particular, it is assumed that a thin cylindrical specimen 



is squeezed between two perfectly smooth, impermeable rigid plates such that 
compressive strain in the axial direction (see. Fig. [l]) is homogeneous. 



Rigid impervious 
smooth plate 



Fit) 
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1 
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Fig. 1 Schematics of the unconfined compression test. 



Let F{t) and w{t) denote, respectively, the variable compressive force and the 
vertical displacement of the upper plate with respect to the lower plate (t is the 
time variable). Then, the compressive strain is defined as 



e{t) 



h ' 



where h is the sample thickness. (Note that compressive strain is negative. 
Let us also introduce the non-dimensional variables 



F : 



t 



HAkt 



(1) 



(2) 



where Ha = As + 2/is, /is and As are the confined compression modulus and Lame 
constants of the elastic solid matrix, respectively, k is the tissue permeability, a is 
the radius of specimen, t is the dimensionless time, F is the dimensionless force. 

Finally, let F{s) and e(s) denote the Laplace transforms with respect to the 
dimensionless time, and s is the Laplace transform parameter. According to |Arm-| 
strong et al. (1984), the following relationship holds true: 



F{s) = se{s)Kis), 



(3) 
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3/0 (Vs) 



Kis) 



8fls h{y/s) 
Ha yfs 



2^is h{y/s) 



(4) 



Ha V~s 

Here, /o and I\ are modified Bessel functions of the first kind. 
By applying the convolution theorem to Eq. ([3|, we obtain 



dr 



K(i - t) dr, 



(5) 



where K{i) = jC~^{K{s)} is the original function of K{s), t = 0— is the time 
moment just preceding the initial moment of loading. In deriving Eq. ([5|, we used 
the formula 

C-'{s-e{s)} = ^+e{0+)5(i), 
at 

where 5{t) is the Dirac function. 

By using the results obtained by Armstrong et al. ( 1984j), we will have 



K{t) = 2(1 + Us) + J2 2(1 + ^s)A„ exp(-a^t), 



(6) 



Here, a„ are the roots of the characteristic equation (with Jo and Ji being Bessel 
functions of the first kind) 

J (1-2^.) Ji(a;) _ 

■JOKX) - -7- r U. 

The inverse relation for Eq. ([5| can be represented as 

t 



dF{T) 

dr 



M (i - t) dr, 



(8) 



where M(t) = C~^{M{s)} with M(s) defined by the formula 



M{s) 



Ha y/s ) 



Again, making use of the results by Armstrong et al. (19841, we get 

1 ^ B„ 



M(t) 



Bn 



2{l + i^s) 



^ 2(1 



\ 2(1 + 1/,) 



exp(-/3^£), 



9(1-^^)2/32 -8(1 + i'.)(l-2i/.)' 



(9) 



(10) 
(11) 
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where are the roots of the corresponding characteristic equation 

A{1 - 2iys) Ji{x) 



Jo{x) 



3(1 - Vs) X 



0. 



The short-time asymptotic approximation for the kernel K{t) (the same can be 
done with M(t)) can be obtained by evaluating the inverse of K{s) as s — ^ oo. For 



this purpose, we apply the well known asymptotic formula (see, e.g., (Gradshteyn 
and Ryzhik||1980[ )) 



1- 82; 



(12) 



Making use of ( 12 ) , we expand the right-hand sides of ^ and ^ in terms of 1/ \/~s. 
As a result, we arrive at the following asymptotic expansions (note a misprint in 
formula (36b) in (Armstrong et al. 1984)): 



Kit) 



M{i) 



2(1 - 2u,) 

V^(l - Vs) 

2(1 - 2v,) 
9v^(l - vs) 



'i+0{i), i- 



(13) 



(14) 



The asymptotic approximations ( |13[ ) and (14) can be used for evaluating im- 
pact unconfined compression test with the impact duration relatively small com- 
pared to the so-called gel diffusion time for the biphasic material tg = a? /[HAk). 

In the dimensional form, Eqs. ([s]) and ([s]) can be recast as follows: 



F{t) 



dw(T) 
dr 



K{t - t) dr. 



(15) 



'At) 



h 



t 

I 



dF{T) 
dr 



M{t - t) dr. 



Here we introduced the notation 

m = 

M(t) = 



Pn 
t 



1 + exp 

n=l 

1- ^B„exp(- — ), 



oiHAk 



(16) 

(17) 

(18) 
(19) 



Note that the sequences pi > p2 > 



a 

and ri > T2 > . . . , which are defined by for- 



mulas (19), represent the discrete relaxation and retardation spectra, respectively. 
Recall that the coefficients An and _B„ are defined by formulas ([T]) and (11). 
Further, in view of (13) and (14), we will have 



Ko = K{G) 



2{i + v,y 



(20) 
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Mo = M(0) = -{l + i^s). 



Hence, the following equalities take place: 



1 + A„ 



2{l + i^,)' 



Now, using formula (21), we can rewrite the function (18) as 



M{t) = Mo + ^ B„ h - exp 



(21) 

(22) 
(23) 

(24) 



By analogy with the viscoelastic model, the functions K{t) and M{t) will be 
called the biphasic relaxation and creep functions for unconfined compression. 



3 Cyclic compressive loading 



Following Li et al. ( 19951), we assume that a biphasic tissue sample is subjected to 



a cyclic displacement input 

w{t) = [wo{l — cosujt) + wi]H{t), 



(25) 



where H{t) is the Heaviside function, wi/h is the prestrain resulting from the 
initial deformation applied to the sample to create the desired preload, wq is the 
displacement amplitude, i.e., wo/h is equal to one-half the peak-to-peak cyclic 
strain input superimposed on the prestrain, and uj = 27r/ is the angular frequency 
with / being the loading frequency. 



Differentiating ( 25 I , we get 

dw{t) _ 
dt ~ 



H{t)wQU! smuit -\- wi5(t). 



(26) 



Now, substituting the expression (26) into Eq. (15), we arrive, after some algebra. 



at the following resulting stress output: 



m 



Es 



wiK(t) 



Wo [Ki (w) COS Lut — K2 {ui) sin ujt\ 



Here we introduced the notation 



2 2 a 

Pn^ An 

+ Pl^'' 



(27) 



(28) 
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PnUjAn 



(29) 



To assign a physical meaning to the introduced functions K\{u}) and K2(ui), 
let us compare the oscillating part of the input strain, that is — {wo/ h) cos u)t, 
with the corresponding oscillating part of the compressive stress, which is equal to 
—Es{'wo/h) [Ki{uj) cos — K2{uj) sinwt] . By analogy with the viscoelastic model, 
we obtain that Ki{lj) and K2{tjj) represent, respectively, the apparent relative 
storage and loss moduli. Correspondingly, the apparent loss angle, 6{oj), can be 
introduced by the formula 

The apparent loss angle S{lo) describes the phase difference between the displace- 
ment input and force output. 



In the case of load-controlled compression, following Suh et al. (1995), we will 



assume that the tissue sample is subjected to a cyclic compressive loading 

F{t) = [Foil - cosiot) + Fi]H{t), 
where Fo is the force amplitude, and Fi is the initial preload. 



(31) 



After substitution of the expression (31) into Eq. (16), we finally obtain the 
following resulting strain output: 



w{t) 



1 



■na^Es 



FiM{t) 



/ °° 2 2 o 

,77/1 '^"'^ 

V 1 1 + T^UJ^ 



Fo[Mi(uj) cos u]t -\- M2 (w) sin wt] 



(32) 



Here we introduced the notation 



°° 2 2 D 



(33) 



00 „ 
M2{u;) = y 



(34) 



Note that Mi{uj) and M2{uj) have a physical meaning of the apparent relative 
storage and loss compliances, respectively. 



4 Displacement-controlled unconfined compression test 

Consider an unconfined compression test with the upper plate displacement spec- 
ified according to the equation 



w{t) = Wo sinijjt, t G (0, tt/cli). 



(35) 
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The maximum displacement, wq, will be achieved at the time moment tm = 
it/{2uj). The moment of time t = t\j, when the contact force F{t) vanishes, de- 
termines the duration of the contact. The contact force itself can be evaluated 



according to Eqs. (15) and (351 as follows 



F{t) 



TTO^Es 



t 

wqU! J COS CUT K{t — t) dr. 



(36) 



According to Eq. (361, the contact forced at the moment of maximum displace- 
ment is given by 



where we introduced the notation 



TTO^Es 



(37) 



7r/(2w) 

Ki{uj) = uj J K{t) sinujT dr. 





(38) 



By analogy with the viscoelastic case, the quantity Ki{uj) will be called the in- 
complete apparent storage modulus. 



Substituting the expression (171 into the right-hand side of Eq. (|38|), we obtain 

UjpnAr. 



+ 1 



ujpn - exp 



(-—))■ 



(39) 



Now, taking into consideration Eqs. (281 and (391, we may conclude that the 



difference between the apparent storage modulus Ki{u}) and the incomplete appar- 
ent storage modulus Ki^oj) is relatively small at low frequencies. To be more pre- 
cise, the difference Ki{lo) — Ki{uj) is positive and of order 0(cjpi exp(— 7r/(2ci;pi))) 
as oj — > 0, where pi is the maximum relaxation time. 

In the high frequency limit, the upper limit of the integral ( |38[ ) tends to zero as 
UJ increases. Thus, the behavior of Ki(ijj) as oj — > -|-oo will depend on the behavior 
of K{t) as t — > 0. According to ( 13 1, as oo, we will have 



(40) 



where Si/2 = Jq ^/xainxdx. 

On the other hand, due to (281 
holds true: lim_fCi((jj) 
Ki{io) ~ Ki{io) for UJ 



20| , and (221, the following limit relation 



Ko as — ^ 00. Thus, in view of (40), we conclude that 
oo as well as Ki(uj) 



Ki(uj) for 0. In other words, 
the incomplete apparent storage modulus Ki (uj) obeys both asymptotic behaviors 
of the apparent storage modulus Ki{uj). 



8 



I. Argatov 



5 Force-controlled unconfined compression test 

Consider now an unconfined compression test witii tiie external force specified 
according to the equation 

F{t) = Fosmcjt, te{0,7r/uj). (41) 

The maximum contact force, Fq, will be achieved at the time moment Im = 
7r/(2w). The moment of time t'j^ = tt/oj, when the contact force F{t) vanishes, 



determines the duration of the compression test. According to Eqs. ( |16[ ) and (41 1, 
the upper plate displacement can be evaluated as follows: 



w{t) = — ^^Foio I cosuTM{t - t) dr. (42) 



Due to Eq. ( 42 I , the displacement at the moment of maximum contact force 
is given by 

w{tM) = — J-^FoMi(c^), (43) 
where we introduced the notation 

Mi(uj)=u: I M{T)smojTdT. (44) 
b 

By analogy with the viscoelastic case, the quantity Mi(lo) will be called the in- 
complete apparent storage compliance. 



Substituting the expression ( 18 1 into the right-hand side of Eq. ( 44 ), we obtain 



M,(.) = Mo + E (l + -p(-£-) ) • (4^) 

71 = 1 



In the same way as it was done with A'i(aj), it can be shown that the incom- 
plete apparent storage compliance Mi(tj) obeys both asymptotic behaviors of the 
apparent storage compliance Mi{uj), that is Mi (a;) ~ Mi(cj) for (jj — > as well as 
Mi(a;) ~ Mi(a;) for lo ^ oo. 
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